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Abstract 

Stock markets can be characterized by fat tails in the volatihty distribution, clus- 
tering of volatilities and slow decay of their time correlations. For an explanation 
models with several mechanisms and consequently many parameters as the Lux- 
Marchesi model have been used. We show that a simple herding model with only 
four parameters leads to a quantitative description of the data. As a new type of 
data we describe the volatility cluster by the waiting time distribution, which can 
be used successfully to distinguish between different models. 



1 Introduction 



Financial stock markets exhibit several universal properties, socalled 'stylized 
facts' [1,2]. There are no time correlations in the index for time scales larger 
than a few minutes [3,4]. The distribution of relative changes or returns are 
power law distributed at large values with an exponent in the order of 3-4 
[4,5]. The absolute returns or volatilities exhibit multifractal moments [6]. In 
the time development of the volatilities we find two features. The correlation 
decays very slowly. Wether this decay is described by a power law with a 
small exponent or an exponential cannot be decided. A power law would be 
equivalent with a 1// noise [4] . The Joseph effect [7] manifests itself by volatil- 
ity clusters, that large changes of the index occur preferably at neighbouring 
times. 

These facts are interesting enough to motivate also physicists to look for mi- 
croscopic models for an at least partial explanation. By microscopic models 
we mean a market consisting of agents following certain trading strategies to- 
gether with a mechanism for the evolution of the price or the index. Many 
different approaches have been persued in the recent years (for reviews see 
[8] or recently [9]). All these models use inhomogenous agents in contrast to 
the effective market hypothesis [3]. The agents are distinguished by differences 
in the behaviour as noise traders or fundamentalists [10-13], the wealth [14], 
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the bidding prices [15], the trading thresholds [16] or the length of memory 
[17]. Alternatively equal agents select different strategies depending on an util- 
ity function [18-21]. The mechanism for creating volatility clusters may be a 
memory, nonlinear coupling between the price and agent parameters or the 
herding effect. The latter may be acchieved by a next neighbour interaction 
as in statistical mechanics [16,19-21] or explicitely in the dynamics [11,12]. A 
model including all three mechanisms is the Lux-Marchesi model [12] (here- 
after abbreviated by LMM). If such a model accounts successfully for the 
styhzed facts as it happens for LMM one would hke to know whether all three 
mechanisms are needed or one alone is already sufficient. In the present paper 
we investigate the question whether herding alone can explain the stylized 
facts. 

As a model we use a simplified version of the LMM. As pointed out in [11] 
there must be three type of agents for which we take noise traders (pessimists 
and optimists) and fundamentalists as in LMM. An important unique feature 
of the LMM is that agents can change their oppinion. Herding is described by 
transition rates proportional to the number of agents in the new state. If we 
neglect the dependence of the agent numbers on the price we can use the ant 
formalism of Kirman [22] or its generalization by Aoki [23] . Since the price has 
disappeared from the dynamic we have to find a phenomenological relation in 
terms of the agent numbers. For that we use the observation [24], that in the 
LMM the dynamics of the price is much faster than that of the agents and 
therefore the equilibrium relation can be used. Assuming a constant number 
of fundamentalists the Kirman model can acchieve qualitative agreement with 
the data [24]. However, we will show that the time structure of the clusters is 
in disagreement with the data. Moreover there exists a winning strategy for 
the fundamentalists, since noise traders are not allowed to become fundamen- 
talists. Therefore we use a three component herding model where herding is 
included also for the fundamentalists as in [11]. 

The paper is organized in the following way. In the next section we describe 
the data (German DAX index [25]). We use the volatility distribution, the the 
correlation function of the volatility and as a new feature the waiting times 
between volatilities larger than a certain minimum value. Section 3 is devoted 
to the description of the agent strategies and the price relation. For large num- 
bers of agents the Kirman model can be written in terms of a Fokker-Planck 
equation. Instead of using transition probabilities we formulate the dynamics 
by the corresponding Langevin equation for the fraction of agents. This model 
and its consequences are described in section 4. The formalism can be easily 
generalized for three components (section 5). Our conclusions are contained 
in section 6. 
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2 German DAX index 



The indices p{t) of stock markets show a long term increase. This is a rather 
small effect, f.e. for the German DAX index 0.03%/day over 30 years. Apart 
from that there are no time correlations on time scales larger than a few 
minutes [3,4]. In contrast the absolute values of the return or the volatility 
defined by 
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exhibits the stylized facts mentioned in section 1, if we consider dayly varia- 
tions. Unfortunately also the values of v{t) averaged over few years show long 
term variations, as pointed out in [26]. To avoid a long term bias we chose the 
dayly quoted values of the DAX in the time period 1996-1999 corresponding 
to 784 trading days, where the three month average of v{t) is compatible with 
a constant value Vav = 0.011. In Fig.l we show the distribution of v/vav on a 
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Fig. 1. The distribution of normalized volatilities for the DAX is shown on a double 
log scale (data points). The solid (dashed) line show the prediction of the three 
component model (Lux-Marchesi model). 

double log scale for these data. The use of normalized volatilities [4] is impor- 
tant when comparing the data with models. Since Vav defines a time scale, the 
translation of computer time steps in a simulation into a real time is taken 
into account using the normalized volatilities v/vav As one sees from Fig.l 
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both LMM and the three component model to be discussed in section 5 are 
compatible with the data. It is difficult to decide whether there is a power law 
behaviour at large v or not. The errors of the data points shown in Fig.l (and 
in the subsequent figures) are computed by binning the volatilities and taking 
•\/]Vb as statistical error, if the bin contains events. 

We quantify the effect of volatility clustering by the waiting times r between 
occurences of volatilities larger than a minimum value Vmin measuring the 
time t = rAt in units of At = 1 day. If the volatilities v{t) are statistically 
independent the probability p^(r) for a waiting time t — 1, 2, ... is given by 

PAr) = -^r-^Y . (2) 

To - 1 V To / 

where the average waiting time Tq is the inverse of the cumulative probability 
P{v > Vmin)- In Fig.2 we show p^(r) for different Vmin- For Vmin « Vav it 
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Fig. 2. The distribution of the waiting time for a volatility larger than Vmin for 
Vmin = 0.004 and 0.016. The curves correspond to independent volatilities and are 
given by equ.(2) using the observed average waiting time tq. 

agrees with (2) indicating that small changes of the index can be considered 
as independent. For values Vmin ~ we see an excess at small and large 
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value of T and a depletion near r ~ tq as expected from formation of volatility 
clusters. The experimental distribution can be characterized by the average 
waiting time To(vmm) and a clustering parameter j^v) defined by 




livmin) = sign(To - t) — Pw{r) ] / , (3) 




where N{t) corresponds to the observed number of events at time r and A^^o to 
the total number. The clustering parameter 7 is proportional to the difference 
of the areas between the data and the curve ]9^(r) expected from independent 
volatilities in the intervalls r < Tq and Tq < r < 2tq. The average waiting time 
Tq{v) is compared in Fig. 3 with independent Gaussean distributed volatilities 
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Fig. 3. The average waiting time tq{v) for volatilities larger than v as function of the 
normalized volatility on a double log scale. The solid (dotted, dashed, dotted-dashed) 
line show the prediction of the three component model (Kirman model, random 
walk, Lux-Marchesi model). 

given by 

- = ^ erfc(^^) (4) 
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and the LMM, which both are remarkably close to the data for v < 2vav 
Again we used the scaled values v/vav of the minimum volatility v = Vmini 
which allows an absolute prediction for one parameter models as in equ. (4). 
In Fig. 4 the cluster parameter ^{v) from equation (3) is shown as function of 



cluster parameter 7 




Fig. 4. The data show the cluster parameter 7 defined by equ. (3) for the DAX as 
function of the normahzed volatility. The solid (dotted, dotted-dashed) line show 
the prediction of the three component model (Kirman model, component model 
(Kirman model, Lux-Marchesi model). 7 = corresponds to no clustering. 

v/vmin- The linear increase of 7(fmm) shows that there exists clustering on all 
scales of v. In the LMM and the Kirman model to be dicussed in section 4 the 
cluster are generated by an intermittency effect by two almost stable states. 
In the second model the fluctuations outside the cluster are almost zero and 
during the cluster of the same size leading to constant 7 and Tq (see Fig. 3). 
The reversed happens in the LMM. Here the fluctuations are independent for 
V < Vav and the clustering begins only at large values of v. Note that the 
LMM has not been adjusted to these data. In view of the many parameters 
in this model the disagreement should not to be taken too seriously. As last 
experimental quantity we will use the correlation function deflned by 
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A meaningful experimental determination of C(r) requires much more statis- 
tics than used in the previous figures. Therefore we use all DAX data from 
1969 to 2002. The correlation function for the volatility is shown in Fig.5 on 
a log-hnear scale. The data show a strong decrease at small values of r which 
compares well with the prediction of the LMM. It fails to account for the 
values at large r, which in turn are in agreement with the three component 
model. 
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Fig. 5. The correlation function for the DAX (crosses) on a log-linear scale. An 
exponential decay would correspond to a straight line. The solid (dashed) line gives 
the prediction of the three component model (Lux-Marchesi model). 



3 Equation of state for the price 

In order to explain the features of the DAX discussed in the previous section 3 
we introduce as in the LMM[12] N agents, which can adopt one of the following 
strategies. Optimistic (pessimistic) noise traders buy (sell) at each time one 
unit of stock. Fundamentalistic traders change their stock proportional to the 
difference of the actual price p(t) and an assumed true value pj. The system 
is characterized by the numbers nj(t) of agents following strategy i = -|-, — / 
(optimistic, pessimistic and fundamentalistic) and the price p(t). In general 
the dynamic will couple Uiii) and p(t) in a complicated way. The new price 



7 



p{t+Atp) will depend on the demand D and the supply S. To avoid dependence 
on splitting orders it should depend only on the difference u — D — S which 
is given in our case by 

Tl f 

ijj{p) = —{Pf - p) + n+ - n_ . (6) 
Po 



The following fairly general form of a price dynamic 

p(i + Ag=pW9(j(^j^) (7) 



has been discussed by Helbing et al.[27]. In relation (7) \{t) can describe a 
reservoir of stocks to fulfill any order or the market liquidity. The montonically 
increasing function g must satisfy g{l) — 1. For the following it is irrelevant 
wether one assumes an exponential form [28], a linear function [12] or a power 
law[27]. Motivated by the simulation of the LMM we assume the price changes 
due to to the dynamic (7) occur at a much smaller time scale Atp than changes 
of ni{t) or \{t) at the scale At. Therefore the dynamic (7) will reach a fixed 
point before any changes in ni{t) or X{t) occur. The mapping (7) has the stable 
fixed point u!{p) — if the stability condition 



i-s'(i)^^ 

Po A 



< 1 (8) 



holds, which can be achieved by sufficient large liquidity A. From the fixed 
point condition a;(p) =0 we obtain the equation of state for the price 



P = Pf+Po ■ (9) 



This relation has been also proposed by Alfarano and Lux [24]. Alternatively 
one can consider (9) as a phenomenological relation expressing that high (low) 
price are expected, if the optimists (pessimists) and a price p/ if fundamental- 
ists dominate the market. In order to avoid negative prices the second term in 
relation (9) must be much smaller than pj. Since the average returns for the 
DAX or other markets are in the order of 1% the probabifity for a negative 
price is very small. Neglecting quadratic terms in v we find for the volatility 

(1) 



vit) 



Po 

Pf^f 



An^ An_ 



H 

At 



At 



(10) 



8 



4 Kirman model 



In this section we simplify the model even further by keeping the number of 
fundamentalists constant. This serves mainly to explain the method. It can be 
easily generalized to three components (see section 5). It has been introduced 
by Kirman [22] for volatility clustering and behaviour of ants. The noise traders 
follow the strategies denoted by i = +, — . The transition rate n per unit time 
for one agent to switch from strategy i to k is given by 

n{i k) = ttik + hik ■ Uk . (11) 



The number of noise traders 



_+n_ is to be kept constant. The constants 



a describe a spontaneous rate and b correspond to the herding behaviour. For 
symmetry between optimistic and pessimistic traders a and b do not depend 
on i, k. The property bik — bki has the remarkable consequence [22,29] that 
the equilibrium distribution wo(?i+) exhibits non thermodynamical behaviour 
in the sense that it becomes for large ric a non trivial function of 



X 



n+ — n_ 



:i2) 



instead a Gaussean around the mean value of x with a width oc 1/ This 
solves a problem in the LMM and other models which become trivial in the 
limit of many agents [30]. The transition probabilities (11) lead to a master 
equation. In the limit of large Uc it can be written in the form of a Fokker- 
Planck equation for the distribution w{x,t): 



Aw(a;, t) 



d_ 
dx 



1 d 

-A^{x) w{x,t) + 2-Q^Dx{x) w{x,t)) 



(13) 



with a drift term 



Aj.{x) — —2ax 



(14) 



and a diffusion term 

D,{x) = 2b{l - x^) 



(15) 



The equilibrium distribution wq{x) depends only on the ratio e 
rate constants 



wo{x) 
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with the normahzation given by Euler's S-function. If a system is described 
by a Fokker-Planck equation, we can replace the original dynamic (11) for the 
individual agents by a Langevin equation for x 



Ax 
At 




(17) 



with a white noise r]f In equations (13) and (17) no dependence on appears. 
It is hidden, however, by the condition, that the time scale and the rates 
must satisfy 

aAi, 6Ai oc l/ric « 1 , (18) 



if the dynamic (17) should represent the original dynamic (11). In the numer- 
ical simulation we use (17) in the form 

x{t + 1) = x{t) - epx + sjpil - a;2) r^^ (19) 

with an effective rate p = 2bAt. Using relation (10) for the volatility, we arrive 
at a model similar to the one proposed by Alfarano and Lux [24] 

v{t) ^Vo\x{t + l) - x{t)\ . (20) 



Our version differs from [24] that in the latter the dynamic (11) is used and 
the drift term is replaced by a reflecting boundary at = 1 as in the LMM. 
For small x the system is described in both cases by a random walk. Near 
|a;| = 1 the diffusion vanishes and the system remains a long time near the 
boundary. The average of the return time Ty^ can be computed analytically 
[31]. For small e we obtain 



Leaving after Tyj the boundary x"^ — 1, volatility clusters arc formed by the 
diffusion mechanism (19) with a time scale 1/b << T^. In Fig. 6 the results 
of a simulation are shown with e = 0.02 and p = 0.02. The left part gives 
x{t) for 4 ■ 10^ time steps. The accumulation of points near |a;| = 1 is clearly 
seen. In the middle part we show the distribution of the volatility as function 
of v/vav The behaviour at large v is compatible with the data, but not the 
peak at small v. On the right we show the correlation function of the volatility 
on log-linear scale. It agrees with the observed correlation (sec Fig. 5), that a 
sharp drop at small r is followed by an exponential decay as indicated by the 
straight line in Fig.6. However, the drop is much to small and the decay too 
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Fig. 6. Simulation of the Kirman model. The left part shows the time serie x{t), the 
middle part the volatility distribution and the right part the correlation function 
given by equ. (5). The straight line corresponds to an exponential fit exp(— 7r) with 
7 = 0.022. 

strong as compared to the data. One could get agreement with e values near 
0.5, which would give totally unrealistic volatility distributions. Apart from 
the peak at small volatilities and the quantitatively wrong correlations the 
Kirman model exhibits further unwanted features. Since the diffusion is fast 
the value of the volatilities in a cluster are always of the same size Vq given by 
(20) and do not reproduce small size clustering seen in Fig. 3 and Fig.4. The 
index values are bounded by 

\P-Pf\<- . (22) 



Such a bound apart being unrealistic leads to a winning strategy for the fun- 
damentalists in contradiction the efficient market hypothesis [3]. These defi- 
ciencies can be cured by using a three component model, where noise traders 
can change into fundamentalists. 



5 Three component model 



The three component model is obtained, if one includes in the agent dynamic 
(11) also the fundamentahsts (i, k = -|-, — , /). The system is described by the 
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fraction x as in equ. (12) and the ratio of noise traders to fundamentalists 
u^'^±±^ . (23) 



The total number A?" = n+ + n_ + n/ is to be kept constant. In order to 
obtain a nontrivial distribution in u in the limit of large N we require the 
a symmetric herding matrix hit- The rate parameters should be symmetric 
between optimistic and pessimistic noise traders leading to 

— b+k , Q+fc — O'-k , flfe+ — flfe- ■ (24) 



There remain only five independent rate constants in the matrices a and b. 
A last condition comes from an integrability condition for the existence of a 
master equation, as discussed by Aoki [23]. This reads in our notation 

^ = . (25) 



Now we perform the same steps from the transition rates to the Langevin 
equation for x and u as in the Kirman model. The equation for x is the same 
as equ. (19). For u we obtain a second equation 



u 



(t + l)^ u{t) + AM + \/Du{u) Vt (26) 



with a drift term 

AM = (1 + k) [(1 - es/2)u + e] ps (27) 

and a diffusion term 

Du{u)=psu{l + u)\ (28) 

where we introduced in analogy to (19) the effective rate p3 = 2b^fAt and the 
ratio €3 — a^f /b+f. The expression for the volatility is given by 

v{t) = vq \u{t + l)x{t + 1) - u{t)x{t) I . (29) 
The equilibrium distribution for the fractions x and u factorizes into a product 

, , (l-X^y-' M^c-l^ ^ ^^)-(63+2.) 
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In the two component model volatility clusters are due to a large return time 
for e — 0. Unfortunately, small e are responsible for the unwanted behaviour 
oc in the volatility distribution. In the three component model one can 
avoid this drawback by a large time for a change from a small value Uq to a 
large value ui given by 

T{uo ^ ui) (X {1 + uiy'-\ (31) 

whereas the reversed time T{ui — > uq) is independent of Ui. For simulations 
of the model the parameters e and can be determined from the equilibrium 
distribution of v. Our choice e = 2 and £3 = 6 give a reasonable description 
of volatihty distribution as can be seen from Fig.l. With the rate parameters 
p — 0.001 and ps — 0.01 the empirical average waiting times shown in Fig. 3 
and the cluster parameter 7 shown in Fig. 4 can be reproduced. Especially 
the u dependent return time (31) leads to clustering also for small v, as seen 
in Fig. 4. With these values the correlation function shown in Fig. 5 can be 
predicted, which agrees with the data except in the transition region between 
the drop and the slow decay. Since we are using normalized volatilities, only 
the ratio p/ps = 0.1 is important for waiting times and volatility distributions. 
It means that changes in the index are mainly caused by transitions between 
noise traders and fundamentalists, and not by changes between optimists and 
pessimists. The ratio ea/e = 3 imphes that changes from noise traders to 
fundamentalists are more influenced by the spontaneous rate than the reversed 
process. Since the spontaneous rate may represent private information, this 
seems not to be unreasonable. From the average value oi < u >= 0.4 we find 
that the fraction of noise traders is only roughly 30%. 



6 Conclusions 

The empirical data for a stock market as the DAX can be explained by a 
model based on herding alone. Especially the time structure of the volatility 
clustering can rule out a two component model. At least three components 
are needed. Since the description of the data is comparable or better than the 
predictions of the Lux-Marchesi model, one may conclude that herding is the 
dominant mechanism in the LMM, and the nonlinear coupling of the price 
evolution to the transition probabilities is less important. 
From the values of herding parameters (pa ~ lOp) we find the surprising fact 
that the herding effect is much stronger for an agent in the fundamentalistic 
state than in the noise trader state. A possible interpretation may be the fol- 
lowing. A typical fundamentalistic trader may be a fond manager. She follows 
the present majority if herding dominates. In the opposite strategy not to fol- 
low the trend a heavy risk for her career is buried. If the majority turns out to 
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be successful she will belong to the few who experience losses and will lose her 
job with no opportunity to revise her strategy. This does not happen if she 
has followed the majority. On the long run agents not following the general 
trend will die out. Therefore there exists a strong tendency towards herding 
behaviour also for the fundamentalists. 
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